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ABSTRACT 



This dissertation reports the interaction of intense radiation 
with iodine vapor, neodymium-doped glass, and benzene. A new stimu- 
lated effect was observed in benzene. Broad-band stimulated emission 
in the region of 910oX to 10,900& was observed when the benzene was 
pumped with light at 6943& from a ruby Q-switched laser. The emis- 
sion consisted of a pulse of less than 10 nsec duration. Peak power 
of the pulse was approximately 500 KW when the ruby laser output was 
40 MW. Threshold was estimated at 25 ± 10 MW/cm 2 . 

Stimulated emission at 1.06 microns was observed from neodymium- 
doped glass when it was pumped with Raman-shifted light at 7458&, 
8053&, and 8751&. The Raman-shifted light was obtained by passing 
light from a Q-switched ruby laser through benzene. The emission 
consisted of a pulse of less than 10 nsec duration. 

Absorption of laser radiation at 6943& and 5310& was measured 
in iodine vapor but no stimulated effects were apparent. 
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I. INTRODUCTION 



Shortly after the first laser was operated successfully by 
1 2 

Maiman , Hellwarth proposed controlling the cavity of the laser 
by Q-switching to produce a very intense pulse of light. The 
intensity of the light produced by Q-switching can be so great 
as to produce a large number of non-linear optical effects in 
matter. 

The first of these effects as identified by Eckhardt and 
3 

Hellwatch, et al, was stimulated Raman emission in nitrobenzene, 

benzene, and several other organic liquids. Stimulated Brillouin 

scattering was first reported in quartz and sapphire by Chiao, 

4 

Townes, and Stoicheff. Stimulated Rayleigh scattering and 
Brillouin scattering in nitrobenzene were reported by Cho, Foltz, 
Rank, and Wiggins. Self-focusing of intense light beams, as pro- 
posed by Askargan,^ was found by Chiao, Townes, and Garmire^ to be 
inextricably involved in the other non-linear processes. These 
effects have since been found in numerous other materials and the 
search continues. 

g 

In 1966, Sorokin and Lankard obtained stimulated fluorescence 

in dyes by passing intense light from a Q-switched ruby laser through 

9 

them. In 1969, Briquet, Jego, and Terneaud observed the generation 
of ultra short pulses at 1.06 microns in neodymium-doped glass when 
it was pumped with light which had been frequency doubled from a 



2 

Q-switched neodymium laser. 

This dissertation reports experimental observations on the inter- 
action of intense radiation with iodine, benzene, and neodymium-doped 
glass. New stimulated effects were observed in the case of benzene 
and neodymium-doped glass. Absorption of the laser radiation was mea- 
sured in iodine vapor, but no stimulated effects were apparent under 
the experimental conditions of this investigation. 



II. RADIATION INTERACTIONS 



A. Emission 



A comprehensive study of radiation interactions with matter 
is obviously a topic too broad to discuss in detail here, and the 
variety of ways in which it can be treated are too numerous to 
treat in this paper. We will, therefore, give a brief treatment 
of the interaction of radiation with matter similar to that used 
by Dirac, ^ as modified by Marcuse.^ 

When considering more than a single isolated system, it is 
often convenient to write the Hamiltonian as 



H = H + H. 

o mt 



( 1 ) 



where H q is the sum of the Hamiltonians of the individual systems, 
and is the interaction Hamiltonian. The Schroedinger equation 

is then 



i* . ( „ o + „ int) | T> 



( 2 ) 



We can transform this by a unitary transformation to the interaction 
representation where the state vector is now 



H 



l'T(t) )= | ¥ > 



( 3 ) 



4 



and the interaction Hamiltonian is 

H 



H 



i -T 2 t -i -j* t 

h; = e h H e 
int int 



In this representation, the Schroedinger equation is 



(4) 



H int l r(t)> 



(5) 



This can be integrated to give 



|r(t)>= |r<o )> + |^/ t H^ nt (t) |r(T))dt 



( 6 ) 



It is not the intent of this paper to discuss the solutions of this 
equation. Suffice it to say that in practice the equation is solved 
approximately by an iteration technique. This involves replacing 
|y^(T))with | Yip)) and solving for |f^(t))for a first order approxi- 
mation 

|T(0)> (7) 






1 + T* f “Int W * 



For the second order approximation, we replace | "^(t) ) with |Y'(t)) 
and solve for |T^(t)) 



1^(0 > = 

2 



1+7 T C Hr (t) dT 

in o int 



+ 7 0 n 1 ) HT (T ) dr dx 

(i 2 ° ° lnt 1 int 2 X 2 



1 ^( 0 )) 



This process can be repeated to as high an order as desired. It 
is then convenient to introduce the scattering operator 



S = 



1 + n i "lot < T > dT 



+ - L rC o' \'nt <\> "lot <t 2 > dT - dT - + 



(ift) 



l 



( 9 ) 



and to write 



5 



K(t)>- s l^(0)> 



( 10 ) 



For a problem of absorption, emission, or scattering, we are 
generally interested in the probability that a system which was 
initially known to be in state |^(0)) at t = 0 (normally an eigen- 
state of H q ) will be found in an eigenstate | ) at time t* | ) 
is an eigenvector of corresponding to the energy eigenvalue, E^. 
This probability is 

p f = i <E f int) > i 2 (id 

which, from (10), can be written 



P = | <E | S |T t0)> | 



( 12 ) 



According to our previous choice, | V "XO) ) = | ) was an eigen- 

state of H also. Hence, 
o 



P f = | < E f | S |E.) 



( 13 ) 



If we now take this to a first order approximation, we have 

H H 

1 . o . o 

S = 1 + / C HT _ (x) dx = 1 + / e 1 ^ H. e dx 

ifi o int ±h o int 



( 14 ) 



and 



P f = 



H 



H 



< E fl E i> + T* {,'<^1^ H i„t E i> dT 



( 15 ) 
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We know, however, that two different eigenvectors are ortho- 
gonal, that is, in our case, (E^|E_^) =0. We also know that when 



H | E ) = E | E ) 
o 1 n n 1 n 



(16) 



then 



f(H )|E > = f (E ) | E > 
o n n n 



(17) 



Hence 



P. = I^T < E r I H - - | E - > 

f 1 ih f 1 int 1 l o 



i(E f - E ) - 

f 1 *dT | 2 



(18) 



and by integration 



±5T <E f l “int l E i> L i(E, - E.) 



F ‘* * f X- 

By means of trigonometric identities, this can be put in to the 

familiar form 



P f ‘ 4 l <E fl “int M* 



sin(E f - E i> - h 



(E f - E t ) 



(19) 



( 20 ) 



At this point, it is convenient to introduce the Dirac delta 
function in the normalized form 

! sin(E - E )-t— 

6(E - E ) = lim TT = 

f 1 _t -> oo E f " E i 

2h 

We will be concerned not only with the actual limiting case, t -* 00 , 



( 21 ) 



but for all cases of sufficiently large t. For these cases, the use 
of the delta function is a good enough approximation for our purposes. 



7 



The probability is then 

sin(E - E ) fj- 

P f ' ^ l< E fl “i„t l E l>! 2 ■' (E 1 -e)) *< E f - E i> 

which by use of the properties of the delta function becomes 



( 22 ) 



p f -^rl< E f l «i„t l E i>l 



6(E f - E i ) 



(23) 



This is the transition probability from | E^ ) to |E^), but if 
we are concerned with the total probability of a transition we must 
sum over all possible transitions. 



P = Z? (24) 

This total probability is normally arrived at by assuming the 
density of final states is great enough to introduce the number of 
states dE^ in the interval dE^, and by replacing the sum with 
an integral. Hence 



P = / P f p f dE f (25) 

and 

p -_r |<p fi H i„ t i E i>i 2 { < E f - E i>f>f dE f <26 > 

which gives 

P=i f l< E fl >< lnt M’p, < 27 > 

If we are interested in the transition probability per unit time. 



this becomes 



w=^=-^ |<E f | I E n . > | 2 P 



int 



8 

(28) 



In the case of an electromagnetic wave interacting with a 
bound electron, the interaction Hamiltonian is given by 

e p ^ -> 0 

H int ’ - m <p ' A) + 2S A 



(29) 



where p is the momentum of the electron and X is the vector poten- 

e 2 

tial of the field. The term A 2 is negligible in all cases 
which we will consider, except Rayleigh scattering. 

It is convenient to take the plane wave representation of A 
expressed in terms of creation and annihilation operators, namely 

2 , -ik *r ik *r 

A = E E T*T ~~ e (fle m +ae m ) 

, m/ 2 eo) V ms ms — 

s=l m 



(30) 



ms 



and to express the initial and final states of the system as pro- 
ducts of the electron state and photon state: 



I E_. ) I ^ ) I n , * n o * • • • • • • ) 



(31) 



E f > = |H' f )| n^, , ... n^+i , ... > 



(32) 



where the are photon occupation numbers of all the radiation modes. 
We can now obtain the transition probability per unit time for emis- 
sion by direct substitution into equation (28). The result is 



w = 



T\e 



m ew^v 



<n w + l)p f | <Y 



-ik *r 

- i e • pe V | Y. ) | 
f 1 va v ' l 1 



(33) 



9 



This can be written in an integral form as 
Tie 2 ^ 2 



w = 



m £0) V 
V 



-V -> 

V ► ~ik * r 

(n + l)p c I ’ (e • W.) e V dVl 2 

va f 1 f va i 1 



(34) 



The density of states is obtained by considering emission into a 

solid angle dQ around k . The final electron level is considered 

v 

well defined, so the density will be determined by the density of 
photon states available. This is given by 

2v ,e A < 

Pc (p. ) TZ df) 

f (2 tt) 3 e o hC 

If we now restrict ourselves to spontaneous emission (n^ a = 0) 

ik • r 

and to dipole radiation e « 1 , we obtain 



(35) 



e 2 ^u) 



w 



v 



spe ,2 2 3 £ 

v 4tt e m e o 

o 



Sf\ ! Y £ * e w • W ± dvP da 



(36) 



For the dipole case only, this is equivalent to the more fami- 
liar form 



03 



w 



1 y va I 2 d n 



s P e 4tt 2 e tic 3 fi 

O 



where the dipole matrix element 



(37) 



va * -> -v 

"fi - e 1 If <%o • r)Y i dV 



(38) 



was obtained by starting with an approximate interaction Hamiltonian 
of the form 



H. 



int 



eE • r 



(39) 



10 



where 



E = - 



3 A 

v 

3 t 



( 40 ) 



If we integrate over all possible orientations of the emitted 
radiation, and take into account the possibility of degenerate 
states, we end up with the total probability for spontaneous emis- 
sion at frequency u) : 



v 



0 ) 



<«_.>= " v ° ‘ E ° In” |» df! 

fc a £l 1 

o 

V0 



spe 



( 41 ) 



The dipole matrix element, is still a function of the 

orientation of the atoms. In a gas, the orientation is random; so 
an average is usually performed over all possible orientations. This 
requires the use of an average dipole matrix |y^| which is indepen- 
dent of the solid angle. If we introduce the average radiation life 
time T and integrate over the solid angle, we then obtain 

to 3 g J%- _ z 

x ^ <w spe > = ire *c 3 " ° l y fjj < 42 > 

spe o 



B. Absorption 

It is obvious that absorption cannot be a spontaneous process. 
There must always be an occupied mode in the initial state before 
absorption is allowed. The derivation of the expression for 



11 



he transition probability for absorption closely parallels that for 
spontaneous emission, except that the density of initial states en- 
ters into the equation instead of the density of final states. They 
are of the same form, except that for absorption there is a factor 
of n^ in place of n + 2 in equation (33). Clearly, n^+ 1 = 0 pro- 
duces no absorption, since this means that there are no atoms in the 
initial state. When the equations are developed along similar lines 
as before, we get 



g 2 s h — i vo 

W = ^ — - ( — ) n hi-. 

abs 4ir 2 e e o f ' 



VO 



VO | 2 

fi 1 






(43) 



where is the average number of photons in each mode. 



C. Multiple Photon Processes 



In addition to single photon absorption and emission, it is pos- 
sible to have multiple photon processes. It is necessary to use the 
second order approximation in calculating the scattering matrix to 
account for two-photon processes. Consequently, the mathematics be- 
comes rather involved. A considerable simplification is achieved by 
considering only two radiation modes. In this case, the vector poten- 
tial of the radiation field becomes 



t - faA )' 5 



V' 
o 



-ik *r ik *r 

1 — i ^ (a e 1 + e 1 ) 



(co ) 2 1 

, -ik *r 



ik *r 



/ T I 

(a e + a e 



(co ) 

2 

The initial and final states for two-photon absorption are 



( 44 ) 



IV = IV I n > I n ) 

1 1 1 2 



( 45 ) 
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|E-> = |V | n -l) | n -1) (46) 

1 t t i '2 

and, as usual, the interaction Hamiltonian is 

H = - - p • A (47) 

m 

The probability of two-photon absorption is now given by 

P f = | <E f | S | E ± ) | 2 (48) 

but the first order does not contribute; so, using the second or- 
der, we have 

(49) 

Using steps similar to those we used to find the first order proba- 
bility, and inserting the operator 



p f = 



(e. 



(ift) 



- / dx / 1 H 7 .. (t ) HT .. (t ) dx | E. ) 
2 o i o int i int 2 2 i 



I = E IE ) <E I (50) 

V V V 

where the |E^)vectors form a complete set of eigenvectors for H q> 
we have 



p f = 



— l { <eJ H. |E )<E | H. |E.) 

(i *) 2 V f lnt V v int 1 



X T 

„ X i(E. - E J- 1 - i (E - E.)- 1 

r t , x f v h V x A i 

/ dx / e e dx ) 

o 1 o 2 



(51) 



This can be integrated to give 



p f - 



- 2i 5 { 



< E f H int V < E V H i„t E i> 
E i - E v 



1(E f - Vt* sin(E i - E f>f* 



E i - E f 



( 52 ) 



e 



} 
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Here we have neglected the quantity 

T 



,* e i(Ep ' ^ 
O 



103T 



dT = 



f 

o 



dT 



(53) 



because it is a rapidly oscillating function of t, and hence, gives a 
negligible contribution to the total integral. The transition proba- 
bility per unit time is, accordingly. 



_ dP d 2 tt 

“ ■ di = i? ; p f N f dE f * — 



5‘- 



<e Jh. Je ) (e |h. |e. > 

N f 1 int 1 V f N V 1 int 1 l 

E. - E 
l V 



Np 



(54) 



where W is the number of atoms and p is their energy density. 

If we now insert the terms for H. and A, we obtain the follow- 

mt 



ing result: 



w = 



7T h 



n n 

L_2_ (£)* Wp 

2e 2 v 2 a) 1 u 2 m 



<y f l e |¥> <fl t-f |0 



B { - 



v 

W. - W + *03 
1 V 1 



<'f'J e -ply ) (U e -p If.) 
f 1 i^'v v 1 2 1 i - 



(55) 



- W + hu 
i V 2 

where the T's are electron wave functions, and the W's are the corres- 

-y -y 
ilc * r 

ponding energy eigenvalues. Here it is assumed that e is constant 

over the range where T gives an appreciable contribution to the inte- 
grand. 

Two-photon absorption can be interpreted as the absorption of one 
photon, which raises the system to a virtual state 4^, together with 
the simultaneous absorption of a second photon, which raises the sys- 



tem to the final state, 

If there is only one radiation mode so that = 03^ , then we get 



w = 



7T h n(n - 



2e 2 V 2 



(r) 4 % 



03 



m 






<H' f | e*p |¥ v ) (Vj e*p |’F i )_ 



W . — W + h 03 

l V 



(56) 
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It is interesting to note that even though the probability 
for two-photon processes is much lower than for single photon 

O 

processes, the probability is proportional to n ; that is, the 
transition probability is proportional to the intensity of the 
light squared , whereas single photon processes are linearly depen- 
dent on the intensity. This is the reason that higher order proces- 
ses come into importance when light from high-powered Q-switched 
lasers is passed through matter. 

By now, it should be clear that probabilities for three-photon 
and higher order absorption can be calculated by using higher orders 
of the scattering matrix. Clearly, multiple-photon emission is also 
a possible process in this regime of high-power radiation fields. 

D. Raman Scattering 

The literature on Raman scattering is amazing from the standpoint 
of ways in which the theory is treated. These range from a quantum 
mechanical second order perturbation theory through a semiclassical 
third order perturbation theory to a strictly classical treatment. 

Each of the different methods have value, so we will briefly examine 
two of them. 

From a classical point of view, we can account for Raman scat- 
tering in the following way. Consider a molecule as consisting of an 
electron cloud surrounding the nuclei. Further, let the incident ra- 
diation on the molecule be represented by a monochromatic wave whose 



IS 



electric vector at the position of the molecule is given by 



E = E sin a) t 
o p 



( 57 ) 



The polarization produced by such a field is 



| P I = a| E 



( 58 ) 



where a is called the polarizability of the molecule. This polari- 
zation is due to the positive nuclei moving one direction under the 
influence of the field E, while the electron cloud is displaced 
slightly the other direction. Under the influence of visible or 
ultraviolet light, the frequency is so great that the heavy nuclei 
cannot respond; hence, the polarization is due almost entirely to 
the movement of the electron cloud. 

If the nuclei are vibrating, it is obvious that this will affect 
the polarizability of the molecule slightly, due to the change of the 
relative positions of the nuclei and electron cloud. To a first ap- 
proximation, we can account for this change by writing 



where >> a and U)^ is the frequency of vibration of the molecule. 
Using this value of the polarizability we get 



a = a + a • sin a) t 
o U 



u 



( 59 ) 



P = a E sin a) t + a E sin a) t sina) t 

o o p i o u p 



( 60 ) 
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Using trigonometric identities, this can be rewritten 



P = a E sinoo t + ^ a cos (a) - oo )t - % a cos (a) + u) )t (61) 

oop i pv i p v 

The oscillating polarization classically gives rise to radiation at 
the frequencies of the polarization which are: oo^, corresponding 

to Rayleigh scattering; (co^ - oo^) , corresponding to Stokes Raman 
scattering; and (co^ + U)^) , corresponding to anti-Stokes Raman scat- 
tering. 

The drawbacks to this simple-minded treatment are several. It 
predicts the Stokes and anti-Stokes lines will be of equal intensity. 
It doesn't predict rotational or electronic Raman scattering, and it 
doesn't predict more than the first order lines, although from a 
classical point of view, one would suspect that all of the "beat" fre- 
quencies of u)p and o)^ would be present. The classical point of view 
really isn't very satisfactory. 

Using much of the same procedure as was used in two-photon ab- 
sorption, we can calculate the probability for Raman scattering by 
using the second order approximation of the scattering operator. For 
the interaction Hamiltonian, we will use 



H. = - - (p * A) 
int m 



and the initial and final states will be 



( 62 ) 



|E. > = |T. > | n >|n ) 

1 i 1 i 1 p ' 1 r ' 



( 63 ) 



|E f > = |y f > | n - 1> |n r + 1> 
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(64) 



where is the number of photons in the primary mode, and n^_ is the 
number in the Raman mode. 

The e and p in the interaction Hamiltonian refer to the charge 
and momentum operator of charged particles in the molecule doing the 
scattering. We will omit the sum over all particles for simplicity, 
and remember that it might be an electron or a nucleus. The scattering 
could give rise to the change of electronic, rotational, or vibrational 
state of the molecule. It most frequently results in a change of vibra- 
tional state. 

Performing the calculations as before, we obtain 



TT tl 



n (n + 1) 

ii rt — — 

w " 2e e V z u> co 
p r p r 



<S>V 



3 { 



<V e r *P IV <V e -p IV 



W. - W + hu 
l V p 



<t,| e -p |f ,) <V | e -p |V.> 
fi p r 1 v V 1 r 1 i 



(65) 



W. - W - fta) 
i v r 

For this process, is the total number of molecules in the initial 
state, and p is the energy difference distribution function for the 
final and initial states. 

This transition probability is actually for the Stokes part of 
Raman scattering, i.e., > W^. We can obtain the anti-Stokes part 

by interchanging the indices p and r, and then, since we want to use 
the same two eigenstates of the molecule, we interchange the indices 
i and f. Now, using the conservation of energy 



W. - W. = ft (<i) - 03) 

f i p r 



( 66 ) 
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we note that the term inside the absolute value signs for anti-Stokes 
scattering is just the complex conjugate of the value for Stokes scat- 
tering. When the absolute value is taken they are the same. 

The total transition probability per unit time, which amounts to 
emission of a photon at u)^ , is the probability for Stokes scattering 
with a transition from to minus the probability for anti-Stokes 
scattering with a transition from to T^. 

The total transition probability per unit time is 



— / j£\ 4 

w ~ 2z e V z 'nr or co 



v{ 



-H- [i.n 

go go Li 1 

p r 



(n +1) - W n (n +1) 
p r f r p 



>] 



ttyl ypiv 



+ 



(T I e *p IT ) (T |e *p|T . ) 
f 1 p^ 1 v /N v l r^ l i j 



(67) 



w. - W + fiu) 

1 v p 

The total power scattered into dq radiation modes in the frequency 
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interval dcor and solid angle dft is 
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The number of radiation modes of one polarization is given by 
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If we now make use of the fact that the energy difference distribution 
function, p, is normalized according to the following condition 



/ p(W f - WJ d(W f - W i )= hT p(rtoo p - ftu> r ) doo r = 1 (70) 

and also p is narrow enough to consider all other factors constant in 
the range where it has appreciable value, we can obtain 



